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i-^ . Abstract 

■ The mass dependence of the M3Y-type effective interactions and the effects of 

tensor correlations are examined. Two-body nuclear matrix elements are obtained 
by the lowest order constrained variational (LOCV) technique with and without 
tensor correlations. We have found that the tensor correlations are important espe- 
cially in the triplet-even (TE) and tensor-even (TNE) channels in order to reproduce 
the G-matrix elements obtained previously. Then M3Y-type potentials for inelastic 
CD ' scattering are obtained by fitting our two body matrix elements to those of a sum 

of Yukawa functions for the mass numbers A = 24, A = 40 and A = 90. 
'. PACS number(s): 21.30. -x, 21.60.-n, 24.10. -i 
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1 Introduction 

Q . 

5 . The study of inelastic reactions from fundamental nucleon-nucleon interaction has been 

a major topic for nearly 30 years now following the pioneering G-matrix work of Bertsch 
and collaborators [1]. In particular their M3Y interaction [1] and its density-dependent 
^ ' versions [2] have shown that one can predict quite unambiguously inelastic scattering from 

' fundamental nucleon nucleon interaction. Because of their semi-microscopic approach, 

there has been a lot of interest in these interactions by many researchers [3,4]. 

Recently [5] we have produced a similarly motivated potential for A = 16 nuclei which 
was based on the lowest order constrained variational (LOCV) approach. We have found 
that our interaction was very similar in the singlet-even, triplet-even, tensor-even and 
spin-orbit odd channels with those of the M3Y interaction but differed significantly in the 
singlet-odd, triplet-odd and spin-orbit even channels. It thus appears that the issue of 
finding precise strengths of interactions in these channels has not yet been settled and still 
requires some further attention. The validity of using the M3Y interaction for inelastic 
scattering has been tested in many applications [6] and we believe that it is possible to 
improve our understanding of inelastic scattering from realistic forces if all the correct in- 
gredients such as finding the precise strengths of the interaction in all angular momentum 
channels and its mass dependence are considered in a systematic way. 
Finding these ingredients may not be very easy and may require going beyond the usual 
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non-relativistic two-body approach of nucleon-nucleon forces to include for example, the 
three-body forces and isobar degrees of freedom. Besides these considerations we must 
check carefully our models since these quantities may also be model-dependent. For ex- 
ample, in the G-matrix approach, it was found that the triplet-even (TE) matrix elements 
used in constructing the inelastic potentials depended sensitively on the starting energies 
[1]. This is a problem of mo del- dependence whereas in our model, these same quantities 
depend sensitively on the strength of the tensor correlations. Without tensor correlations 
we shall show in this paper that we have no agreement in these channels with the G-matrix 
approach. 

Another important consideration is the mass dependence. It is useful to consider the 
mass dependence of such an interaction in a systematic way while using it on all regions 
of nuclei. For example, while the M3Y interaction has been very successful in explaining 
inelastic scattering from fundamental forces, it was constructed for the A — 16 nuclei. One 
is not sure how it varies with the mass number, A. This last point on mass dependence has 
been emphasized in the classic work of Wildenthal [7] who found that it was impossible 
to find a mass - independent two-body effective interaction that could explain all sd -shell 
nuclei. 

The purposes of the present work are two-fold: (i) to extend the work of ref. [5] which 
was designed for the A — 16 system by including the mass dependence of our potential 
model for inelastic scattering and (ii) to study the effects of tensor correlations which 
affects particularly the TE and TNE channels. Other ingredients such as the inclusion of 
density dependence, addition of three-body forces and isobar degrees of freedom etc could 
also be considered but we shall leave them out from this paper and consider our approach 
a first approximation. 

The paper is organized as follows: In section 2 we give a brief summary of the method 
used. In section 3, we define the A^A^ interaction to be constructed. In section 4 we 
discuss our findings in the present work. The final section is devoted to the conclusion of 
the paper. 



2 Formalism 

In this section we briefly review the LOCV prescription discussed by Irvine and collabora- 
tors [8,9] for obtaining an effective two-body interaction which is adopted for the present 
work. In this approach the approximation to a non-relativistic nucleon fluid interacting 
through a two-body potential is first approximated by the Hamiltonian 

^o = E-|^v.^+E% (1) 

i i>j 

where Vij is the two-body potential. In the LOGV approach the translationally invariant 
component of the trial wavefunction is defined as 

= UF^, (2) 
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where U is a unitary operator which transforms the system to the centre of mass rest- 
frame, leaving us with only intrinsic quantities such that we do not have to worry any 
further about the spurious centre of mass motion. In eq. (2), F is a symmetric product 
of two-body correlation functions defined as [9] : 

F = I[M^J), (3) 

i>j 

designed to accommodate the effect of the strong repulsion of the nucleon-nucleon inter- 
action, and $ is a multi-dimensional product of two-body wave functions. The explicit 
form of f^i'ij) is given in eqs. (10) -(13) below. 

Our task is to evaluate the expectation value of the Hamiltonian given by 





^0 






T 





(4) 



with the desired accuracy. Due to its multi-dimensional nature, eq. (4) is difficult to 

calculate and we must make approximations. The cluster expansion technique is the 
approximating technique wc adopt. Here we divide the system into clusters. We evaluate 
the energy of each cluster starting with the two-body clusters and then sum over all the 
clusters to obtain the total energy: 

E ^ E2 + E3 + E^ + ... . (5) 

To lowest order we minimize the two-body energy term E2 with respect to the functional 

variations of the two-body correlation functions such that only two-body cluster terms are 
important, whereas contributions of E3 and higher order chistcrs are made small. This is 
achieved if the convergence parameter k, defined by the equation below is small [10]; 

^ = ^ /((I - hmyidn, « 1 , (6) 

where Tq is the mean nearest neighbour distance, defined from the nuclear radius it! such 
that: 

R = r-oA'/', (7) 

while < ... > signifies channel average, inchiding the tensor operator. This approximation 
impUes that we have accepted to work only with two-body clusters. Thus our E2 takes 
the form [10]: 

E2 = I E(/2(U)(4/M + V,j)f2m I (8) 

i>j 

where Pij = -^{Pi ~ Pj) is the relative momentum of the two-particle system; M ~ m]\fA 
is the total mass of the nucleus and Vij is taken to be the Reid [11] soft - core potential. 
Since the Reid [11] potential has the form: 

y^j-Evi^ (9) 

A 
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where in different reaction clianneis A, we liave tlie centrai, spin-orbit and tensor compo- 
nents. For eacfi cfiannel we introduce two-body correlation functions, 

/2(^J) = E/^ (10) 

A 

wliere 

/S = fciu,) + fUn,)L ■ S + /^(r,,)5,, , (11) 
with the tensor operator 

Sij = 3ai ■ fij aj ■ fij - Bi ■ Bj . (12) 

In studies regarding nuclear matter and finite nuclei [9], the two-body correlation 
functions were found to have three features which included the 'wound' induced in the two- 
body wave function by the repulsive core of the N-N interaction, the tensor correlations 
especially in the — ^ Di channel and the meson exchange correction. It was found that 
the most important feature of these was the tensor correlations. Therefore the two-body 
correlation functions of eq. (10) have been parametrized in the form [10,12,13]: 

rij>r, : = (1 - e-^(^^^-^=)')(l + aV)'^,,) , (13) 

where = 0.25/m and f3 = 25/m^^. The parameter, a^{A) represents the strength of 
the tensor correlations and is non-zero only in the - ^Di and - channels. 

It should be noted that this form of correlation functions is clearly an over simpli- 
fication. In reality correlation functions are highly density-dependent and are different 
in different tensor channels. Moreover, they are usually obtained from nuclear matter 
calculations. We have retained this form because of its previous success [9,10,12,13] and 
because we wish to study the effects of tensor correlations on our model calculations where 
we have devised a method of switching on and off the tensor correlations. 

In order to evaluate matrix elements we work in the harmonic oscillator basis. For our 
purpose we have picked out only the effective potential energy terms from eq. (8) i.e., 

E', = {^\^Mtm,Mzj)\^) (14) 

i>j 

Here the only parameters appearing in our calculations are; the oscillator size parameter 
and the strength of the tensor correlations, a^{A). We shall later present our results 
with various combinations of these parameters with and without tensor correlations. Fur- 
thermore following the procedure of Bcrstch and collaborators [1] we have separated the 
relative potential two-body matrix elements into those of various channels. These are the 
singlet-even (SE) and singlet-odd (SO) channels denoted by ^5*0 and ^Pi respectively. The 
triplet-even (TE) and the tensor-even (TNE) components were picked from the coupled 
^5"! - ^Di channels. We next define the triplet-odd (TO), tensor-odd (TNO) and the two 
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components of the spin-orbit force ignoring the quadratic term as [14]: 

V(TO) = V(^Po) + 2V{LSO) + AV{TNO) 
V(TNO) = -^[2V{^Po) - SVi^Pi) + P2)] 

V{LSO) = -^[2VCPo) + SVCPi) - 5VCP2)] 

V{LSE) = \ [V{TE) - 2V{TNE) - V{^Di)\ . (15) 

3 The A^A^ Effective Interaction 

In this section we define an effective interaction which is suited to calculations of inelastic 
scattering. We follow the procedure described by Bertsch and collaborators [1] in which 
the two-body matrix elements are fitted to those of the sum of Yukawa functions with 
different ranges. The potential is divided into the central (c), spin-orbit (Is) and tensor 
(t) components as follows: 

^-{r-ij/Rk) 

Vis = Y,^k-^ — T^-T^-'S', 

k inj/Rk) 

p-ijij/Rk) 

^' ^ V"^^''^'-" '''' 

where the dk are the strengths of the interaction which are determined by fitting the 
oscillator matrix elements of cq. (16) to our two-body matrix elements of eq. (14). The 
ranges k < 4 are chosen to be 0.25, 0.4, 0.7 and 1.414 fm, which are motivated by the 
one-boson exchanges. Apparently the longest range 1.414 fm corresponds to the one-pion 
exchange, while the shorter ranges correspond to heavier mesons such as a, p and uj 
mesons. 



4 Results 

The objective of our present paper is to extend the work of ref. [5] in which we defined an 
effective interaction for inelastic scattering for the A = IQ system. In the present paper 
we incorporate the mass dependence as well as investigate the effects of tensor correlations 
on our earlier calculations [5]. 

In Table 1 we present the relative two - body matrix elements of the effective interaction 
for the A = 24 system for the SE, SO, TE, TO, TNE, TNO, ISO and LSE channels. The 
matrix elements were calculated in a harmonic oscillator basis with fiuj = 45/1"^/^ — 
25A~^/^MeV so that fiu; takes on the value of 12.7 MeV in this case. Based on our 
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observations in [8,12] for the sd shell and earlier calculations for the p shell [10], it was 
found that a^{A) is a monotonic decreasing function of the mass number, A and ranges 
from 0.09 at the beginning of the p shell to about 0.07 for the upper end of the sd shell. 
When extrapolated to the A = 90 region, it takes the value of about 0.06. Using this 
information, the calculation is first performed with tensor correlations switched on, i.e. by 
setting a^{A) = 0.075 appropriate for the A = 24 system. However, when a^{A) = 0.0, 
we have only central correlations. We should observe that the tensor correlations are 
effective only in the coupled channels such as ^5"! - ^Di and - ^^2, while the SE and 
SO channels are not subject to the tensor correlations. 

In Table 2 we repeat the same procedure as in Table 1 for the A = 40 system but 
this time we have used a^{A) = 0.070 which we found most appropriate for this system 
with fiu; — ll.OMeV while for the A — 90 system presented in Table 3 we have used 
a^{A) — 0.060 and hu; — S.SMeV based on our observations above. 

4.1 Central Even Channels 

Prom Table 1, we notice an interesting result for the TE channel when the strength of the 
tensor correlations a^{A) is switched off. With a^(A) turned off the TE matrix elements 
are large and positive in our calculation. However, with a^{A) turned on they become 
large and negative. For an S'-wave ("^Si) channel, this attraction is caused by the second 
order (two pion exchange) contribution with the D-wave {^Di) mixing. The role of the 
pion in this channel has long been appreciated. Recently, it is emphasized even more in 
exact calculations of few body systems [15] and in the relativistic mean field theory [16]. 
The TO, TNO and LSO channels are also affected by the tensor correlations through the 
state. The LSE channel is affected by the tensor correlations through the ^Di state, 
however, these effects are not as dramatic as those observed in the TE channel case. The 
SE and SO channels are not affected by the tensor correlations as indicated in the table. 
For the A — 24 system the G-matrix data was unavailable with which to compare our 
results. 

In Table 2 we notice the similar results for the TE channel when the strength of 
the tensor correlations a^{A) is switched on or off. Here also our calculated TE matrix 
elements are large and positive whereas with a^{A) turned on they are large and negative. 
These should be compared with the G-matrix results of Hosaka et al. [17] for the A — 40 
system. Clearly, with ^^(A) switched on to take the value of 0.075, we obtained a near 
perfect agreement for our calculated TE matrix elements with the G-matrix result of 
Hosaka et al.[17] where the difference between the two calculational procedures for each 
(n',n) values is less than 12 percent. Also as shown in Table 2 except for the (n',n) = (2,2) 
quantum numbers, our calculated SE matrix elements are found to agree with the G- 
matrix results of Hosaka et al.[17] to within 25 percent. 

In Table 3 for the A = 90 system we notice a similar trend as in Tables 1 and 2. The 
same observation for the TE channel is inferred in Table 3. Here we have used small node 
quantum numbers to compare our results with the G-matrix results. We notice again 
that the TE comparison with the G-matrix results is impressive when a^{A) is switched 



6 



on. This analysis supports the importance of tensor correlations in effective interaction 
theories that use correlated basis functions. Therefore in subsequent analysis we shall 
only be concerned with matrix elements when a^(A) is non-zero. 

4.2 Central Odd Channels 

As can be seen from Table 2 for the A — AO system, our calculated SO matrix elements are 
similar to the G-matrix for the lower node quantum numbers but they become substan- 
tially different for higher node quantum numbers where for example, for the {n' , n) = {2,2), 
they differ by over 80 percent. This difference between our calculated matrix elements 
and the G-matrix results in the odd angular momentum channels persists and can be 
clearly seen in the TO channel where our calculated matrix elements are positive while 
the G-matrix results are small and negative. We observe a similar trend to Table 2 in 
Table 3 for the A = 90 system. 

4.3 Tensor Channels 

From Table 2 we see that the TNE channels are best represented. In these channel, our 
calculated matrix elements do not differ from the G-matrix calculations by more than 
10 percent except for the (n',n) = (0,0) and (2,2) channels. The TNO matrix elements 
are also very similar to the G-matrix results with our calculated values slightly smaller 
than the G-matrix values. In Table 3 we notice a good agreement of our calculated TNE 
matrix elements with their G-matrix counterparts for the small node quantum numbers. 
Our calculated TNO matrix elements are also seen to be in good agreement with the 
G-matrix results. 

4.4 Spin-Orbit Channels 

In these channels, we notice as indicated in Tables 2 and 3 that our calculated LSO 
matrix elements are in close agreement with the G-matrix values. However, the LSE 
matrix elements are different from their G-matrix counterparts, although absolute values 
of the matrix elements are small. Here, the G-matrix results are negative while our 
calculated matrix elements are all positive. 

It could be seen that in all channels, the matrix elements are decreasing in absolute values 
with increasing mass number, A. 

4.5 Strength of Interactions 

We have performed a least squares fit of our calculated two - body matrix elements in 
the various channels to those of a sum of Yukawa potentials for the central, tensor and 
spin-orbit forms as given in eq. (16). These forms are consistent with potentials for 
inelastic scattering. The selected ranges of 0.25, 0.4, 0.7 and 1.414 fm were theoretically 
motivated and were chosen so as to ensure the OPEP tails in the relevant channels as 
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well as the short-range part which account for the exchange processes. Our interaction 
strengths are presented in Tables 4, 5 and 6 ior A — 24, 40 and 90 systems respectively. 
It should be noted that the oscillator matrix elements of the radial components of eq. 
(16) decrease with increasing A, i.e., with decreasing huj. This is because wave functions 
extend more as A is increased, while the range of the interaction is finite and fixed. Our 
observation here is that the major part of the mass dependence of matrix elements is due 
to the change in the wave function. In order to further study mass dependence of the 
interaction strength, we need more detailed information of the effective interaction. 

5 Conclusions 

The two calculational procedures presented here raise several interesting observations: 
The first observation is that although there is some model dependence in the calculation 
of the matrix elements in the TE channels (such as sensitively on the starting energies 
used in the calculation [1]), reasonable theoretical methods predict attractive matrix ele- 
ments through a significant contribution from the tensor correlations. This indicates the 
importance of the tensor force in the nucleon-nucleon interaction dominated by the one 
pion exchange. 

Secondly, we notice that while the matrix elements calculated in our present method 
are in excellent agreement with the G-matrix results of, for instance, Hosaka et al. [17], in 
most angular momentum channels, there is a problem in fixing the interaction strengths 
especially in the odd angular momentum channels. This problem was also found in ref. [5] 
for the A = 16 system. This can be seen in the SO, TO and the LSE channels. In fact, 
Bertsch and collaborators [1] had to use the set of Elliott [18] matrix elements in the odd 
angular momentum channels to get a good fit for their M3Y interaction. Hosaka et al.[17] 
have traced this problem to the fact that because the harmonic oscillator wave function 
vanishes near the origin for the odd forces, their exist an ambiguity in fitting the G-matrix 
elements in the odd angular momentum channels. 

The calculated matrix elements show an overall decrease in magnitude as we go from 
^4 = 24 to = 90 nuclei through the 74-dependence of hu; and with the moderate mass 
dependence of the strength of the tensor correlations a'^{A). We have calculated the 
strengths of our interaction for A = 24, 40 and A = 90 nuclei as reported here. Since 
these depend on the mass number A, It will be interesting to see if this mass dependence 
may affect calculations based on this interaction. 

Thus, apart from the necessary inclusion of the density dependence and other effects 
in effective two-body interactions, one should not ignore perhaps the mass dependence 
and the model dependence in such calculations in a systematic way. 
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Table 1: Calculated relative matrix elements for A = 24 using Hlo = 12.7MeV. The first 
entry for a column is the result of the present calculation with a = 0.075. The second entry 
in parentheses is the result of the present calculation with a = 0.0. There were no G-matrix 
calculations for comparison for this system. Notice that the SE and SO channels are not affected 
by the tensor correlations. 
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Table 2: Calculated relative matrix elements for A = 40 using huj = ll.OMeV. The first entry 
for a column is the result of the present calculation with a = 0.070. The second entry is the 
G-matrix calculation with the Paris potential. The third entry in parentheses is the result of 
the present calculation with a = 0.0. 
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2 


1 OA 

— 1.30 


r 

— Z.58 


OAT 

—3.95 


2 






A /I A 

0.40 




-1.1172 


-2.3980 


-3.6559 








0.7073 




(-1.55) 


(-3.13) 


(-4.83) 








(0.49) 


LSE 








LSO 








(D/D) 


n = 


1 


2 


(P/P) 


n = 


1 


2 


n' = 


0.15 


0.19 


0.18 


n' = 


-0.15 


-0.35 


-0.53 




-0.0059 


-0.0090 


-0.0091 




-0.3261 


-0.4751 


-0.5801 




(0.03) 


(0.05) 


(0.06) 




(-0.44) 


(-0.68) 


(-0.86) 


1 




0.25 


0.27 


1 




-0.59 


-0.82 






-0.0278 


-0.02831 


1 




-0.6969 


-0.8546 






(0.08) 


(0.09) 






-1.02 


-1.28 


2 






0.33 


2 






-1.07 








-0.0553 








-1.0522 








(0.12) 








(-1.59) 



Table 3: Calculated relative matrix elements for A = 90 using Hlo = 8.8MeV. The first entry 
for a column is the result of the present calculation with a = 0.060. The second entry is the 
G-matrix calculation with the Paris potential. The third entry with parentheses is the result of 
the present calculation with a = 0.0. 



SE 








TE 








(S/S) 


n = 


1 


2 


(S/S) 


n = 


1 


2 


n — 


—4.58 


—4.51 


A rvrv 

—4.09 


n — 


—5.43 


—4.69 


—3.68 




-4.1471 


-3.8435 


-3.2648 




-6.4335 


-6.2354 


-5.5745 












(1.89) 


(2.82) 


(3.57) 


1 




-4.69 


—4.37 


1 




A AO 

—4.48 


—3.65 






-3.7805 


-3.2868 






-6.3463 


-5.7768 














(3.80) 


(4.63) 


2 






A 1 O 

—4.18 


2 






O A A 

—3.09 








-2.9243 








(-5.3810) 
















(5.44) 


SO 








TO 








(P/P) 


n = 


1 


2 


(P/P) 


n = 


1 


2 


n —0 


0.77 


1.13 


1.49 


/ A 

n —0 


0.27 


A O O 

0.33 


A or" 

0.35 




1.2020 


1.2299 


1.1562 




-0.0256 


-0.0591 


-0.0745 












(0.04) 


(0.06) 


(0.08) 


1 




1.82 


2.33 


1 
1 




0.46 


0.53 






1.5768 


1.6314 






-0.0751 


-0.0843 














(0.14) 


(0.18) 


2 






2.99 


2 






0.66 








1.8566 








-0.0813 
















(0.26) 


iJNb 








INO 








(S/D) 


// = 


1 


2 


(P/P) 


// = 


1 


2 


n — 


—2.41 


—3.58 


—4.40 


/ A 


A A A 

0.29 


A OA 

0.30 


0.27 




-2.7015 


-3.9266 


-4.7104 




0.3756 


0.3895 


0.3591 




(-2.84) 


(-4.20) 


(-5.16) 




(0.32) 


(0.33) 


(0.30) 


1 


—1.52 


—2.78 


—3.83 


1 




0.36 


A o r 

0.35 




-1.6481 


-2.8635 


-3.9028 






0.4831 


0.4875 




(-1.78) 


(-3.25) 


(-4.49) 






(0.41) 


(0.40) 


2 


— 1.01 


— 1.97 


-2.96 


2 






0.36 




-0.9310 


-1.9225 


-2.8453 








0.5364 




(-1.17) 


(-2.30) 


(-3.49) 








(0.42) 


LSE 








LSO 








(D/D) 


n = 


1 


2 


(P/P) 


n = 


1 


2 


n' = 


0.085 


0.102 


0.105 


n' = 


-0.084 


-0.207 


-0.321 




-0.0024 


-0.0069 


-0.0094 




-0.1928 


-0.2860 


-0.3548 




(0.017) 


(0.026) 


(0.034) 




(-0.260) 


(-0.400) 


(-0.520) 


1 




0.145 


0.163 


1 




-0.359 


-0.505 






-0.0148 


-0.0183 


12 




-0.4254 


-0.5296 






(0.042) 


(0.054) 






(-0.610) 


(-0.780) 


2 






0.198 


2 






-0.669 








-0.0352 








-0.6610 








(0.070) 








(-0.980) 



Table 4: Best-fit interaction strength (MeV)for A = 24,;iw=12.7 MeV, Q!=0.075 





Channel 


Ri = 0.25 fm 


i?2 = 0.40/m 


i?3 = 0.70/m 


i?4 = 1.414/m 


1 


SE 


10790 


-3538 




-10.463 


2 


TE 


23209 


-7211 




-10.463 


3 


SO 


-321.4 


2593 




31.389 


4 


TO 


421.4 


807.1 




3.488 


5 


TNE 




-1441 


-4.8 




6 


TNO 




-25.9 


23.7 




7 


LSE 


-3239 


1973 






8 


LSO 


-66.7 


-816.7 







Table 5: Best-fit interaction strengtfi (MeV)for A = 40, ^1^=11 MeV, a=0.070 





Channel 


Ri = 0.25/m 


i?2 = 0.40/m 


i?3 = 0.70/m 


i?4 = 1.414/m 


1 


SE 


10427 


-3414 




-10.463 


2 


TE 


22035 


-6741 




-10.463 


3 


SO 


9376 


969 




31.389 


4 


TO 


-138.9 


798.4 




3.488 


5 


TNE 




-1260 


-19.2 




6 


TNO 




11.34 


21.70 




7 


LSE 


-48365 


5188 






8 


LSO 


-1721 


-555 







Table 6: Best-fit interaction strength (MeV)for A = 90,na;=8.8 MeV, a=0.060 





Channel 


Ri = 0.25/m 


i?2 = 0.40/m 


i?3 = 0.70/m 


i?4 = 1.414/m 


1 


SE 


8452 


-2871 




-10.463 


2 


TE 


17554 


-5311 




-10.463 


3 


SO 


3611 


490 




31.389 


4 


TO 


1402 


183 




3.488 


5 


TNE 




-1204 


-25.2 




6 


TNO 




61.02 


19.02 




7 


LSE 


-40737 


4476 






8 


LSO 


-19239 


1731 







13 



